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Abstract 

> : 

We review some recent results on heavy quark production in high energy hadronic colli- 
'. sions. We will discuss in particular the status of production cross sections for bottom quarks 

' and charmonium states and will present some studies on the production of bottom and charm 

CO . jets, at the inclusive level and in association with electroweak gauge bosons. 

> ' 

oo : 1 Introduction 

i> ■ 
on ■ 

i—{ ' Heavy quark production in high energy hadronic collisions consitutes a benchmark process for 
the study of perturbative QCD [0. The comparison of experimental data with the predictions of 
. QCD provides a necessary check that the ingredients entering the evaluation of hadronic processes 
(partonic distribution functions and higher order corrections) are under control and can be used 
to evaluate the rates for more exotic phenomena or to extrapolate the calculations to even higher 
energies. The estimates of production rates for the elusive top quark rely on the understanding of 
heavy quark production properties within QCD. 

Much progress has been made in the field since the last pp Collider Workshop in 1989, both 
experimentally and theoretically. Being impossible to fit in the time allowed a complete review of 
these developments, I will shortly summarise them here before proceeding with the recent material, 
for the benefit of the unaware reader. 

At the time of the last Workshop, one of the outstanding issues in the field was the disagreement 
between the measurement of the CERN experiment UAl of the B cross section at large pt @ and 
the first theoretical calculation of this distribution at the next-to-leading order in QCD (a.k.a. 
NDE, j3|). No results from the Collider Detector at Fermilab (CDF) were available as yet. Since 
then, new data from UAl have come into nice agreement with the theoretical QCD expectations 
and an independent QCD calculation has been performed by a different group ||, providing 
an important confirmation of the NDE results. 

The dust on the UAl new analysis had barely settled, when the success of QCD in describing 
the B production cross section was put into serious embarassement by the first results from CDF 
H pj|, indicating a significant discrepancy between the data and the same theoretical prediction 
working so well at the lower CERN energies. 
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1 Invited talk presented at the 9th Topical Workshop on Proton-Antiproton Collider Physics, 18-22 October 
1993, Tsukuba, Japan. 
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This disagreement at higher energy, whose possibility had already been anticipated in the 
original NDE article, was justified qualitatively by the inadequacy of a fixed order calculation in 
perturbation theory (PT) applied to processes probing initial state gluons at very small Bjorken-x. 
In this regime, in fact, large corrections of order [a s log(l/x)] n ~ [a s log(s/m 2 )] n appear at higher 
orders in PT. Some of the most interesting theoretical papers in this period have been devoted 
to the development of techniques designed to resum these corrections [§. While it became clear 
that indeed these effects can increase the total B cross section by significant factors at asymptotic 
energies, the increase from 0.63 to 1.8 TeV between the CERN and FNAL energies was seen to 
be insufficient to explain the factor of 3 or more observed rate discrepancy. 

Alternative explanations have been attempted, most notably the possibility || that the gluon 
distribution function as determined from standard fits of low energy data is not known sufficiently 
well to allow a reliable extrapolation to the small values of x involved in B production at the 



Tevatron. While there is here some room for improvement, as shown in ||, it is my view [ 10 1 that 
with the most recent data and analyses of structure functions JTl], |l2j it would be hard to account 
for a change in the expected B cross section as drastic as required by the first CDF data. 

In the attempt to further probe the dynamics of heavy quark production, a full NLO calculation 
of the QQ correlations has also been carried out in the recent past [[L3|]. Quantities predicted by 
this calculation (a.k.a. MNR) are, among others, the p t distribution of the quark pair system 
(pf ) and the shape of the A0 distribution, A0 being the difference in azimuth between the 
quarks in the plane transverse to the beam axis. Measurements of this distribution have been 



carried out both by UA1 jljj and by CDF Ul5| , indicating good agreement with the expectations 
of NLO QCD. In this area, a recent study has improved the theoretical calculation of the mass 
and pf® correlations by including the effects of multiple soft gluon emission from the initial state 
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The corrections to the inclusive pt distribution due to multiple soft gluon emission have also 
been calculated, in an important recent paper by Laenen et al. [[T7|]. Since these effects contribute 
terms of order [a s \og(m/pf )} n , they turn out not to affect significantly the bottom cross section, 
while they can be of importance in the case of the more massive top quark, which is produced 
much closer to threshold and for which soft gluon exchanges are more important. A complete 



discussion of the results for the top can be found in ref . [|18 



Parallel to the theoretical developments, signficant improvements have taken place on the ex- 
perimental side. We can now benefit from the presence of the new detector, DO, whose preliminary 
B results have been shown also at this meeting ||20|| . Furthermore, the upgrade of the CDF de- 
tector with a silicon microvertex detector, capable of reconstructing the secondary vertices from 
B decays, is increasing CDF's ability to identify B's with reduced backgrounds This is fun- 
damental, for example, in separating the component of the J/^'s produced directly from those 
produced in the decay of B states. 

The new preliminary results on the pt distribution of fully reconstructed B mesons, presented 
here by D. Crane |fij, suggest a decrease of the measured cross section in the small p t region 
when compared to what previosuly observed || . I collect in Fig. [I] the available CDF data and 
compare them with the NLO QCD calculation || [5]. I pushed to the limit the input parameters of 
the theoretical calculation, namely the renormalization/factorization scale (set to /i = /io/4, with 
/i = yjm 2 +p t 2 ), and the value of A QCD (set to A4^ oop =275 MeV for the MRSDO PDF set 0). 
This choice of parameters still provides an acceptable upper limit to the theoretical expectation. 
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Figure 1: Integrated b quark p t distribution at 1.8 
TeV. CDF data versus NLO QCD. 



Figure 2: Diagrams for the production of a J/lp- di- 
rect (top) and via gluon fragmentation (bottom). 
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As clear from the plot, if the new data were confirmed in the final analyses by CDF and DO, 
the discrepancy with theory would have disappeared, though at the expense of selecting extreme 
values for the theoretical input parameters. 

The situation with the B cross section at the Tevatron is therefore what I would define as 
fluid: new data from HERA in the small x and large Q 2 region will soon provide more reliable 
parametrisations of the gluon densities, and new data and analysis techniques from CDF and DO 
will further reduce the experimental uncertainties of the measurements. 

In the rest of this report I will concentrate on two items: 

• recent progress in the calculation of charmonium production 

• studies of charm and bottom jet production, both at the inclusive level and in association 
with electroweak gauge bosons. 

The interest in these two subjects arises from different observations. The driving force for the 
former is directly related to the measurement of the B cross sections: B decays are in fact an 
important component of the observed rates of charmonium states, J/if>, x an d ip(2S). In the case 
of the ip(2S), for example, it is currently believed from theoretical calculations pi] , |10|| that almost 
100% of the observed events result from B decays. Large fractions are also expected for the J/ip. 
Independent measurements of the B cross section, and knowledge of the decay branching ratios, 
will therefore provide by subtraction a measurement of the direct sources of charmonium, which 
can then be tested against our QCD production models. 

Interest in the latter subject comes from several sources. To start with, the study of heavy 
quark jets (namely jets originating from the evolution of a promptly produced Q or jets con- 
taining a heavy quark produced during the gluon shower evolution) provides yet another test of 
our understanding of QCD. Furthermore, b-jets produced in association with W bosons provide 
powerful signatures for the detection of the top quark [^2| . It is therefore important to understand 
all possible sources of backgrounds induced by direct production of heavy quark jets. The use 
of different data samples which are known not to be contaminated by a possible top, but which 
do contain at some level bottom or charm jets, allows to explore the performance of the heavy 
quark tagging techniques, and to determine the reliability of the theoretical expectations. Finally, 
associated production of E WK bosons and heavy quarks has been suggested as a possible tool to 
study the heavy flavour content of the proton [^3| . 



A more thorough discussion of all of these questions will hopefully appear soon p4|j . 

Before concluding this Introduction, let me just mention an additional couple of very recent 
theoretical developments in the theory of heavy quark production in hadronic collisions, items 
which I will have no time to cover but which I want to point out to the interested reader. The 
first is the study of the interplay between a finite top decay width and the spectrum of soft gluon 



radiation emitted by the top itself and its decay products f25fl . The predicted correlations are very 



difficult to observe, but perhaps future large statistics of detected top and better confidence with 
new experimental techniques will open new exciting possibilities. The second subject is a recent 
calculation of the resummed NLO high p t behaviour of the b cross section, performed by 
matching the NLO fragmentation function formalism |l| with the NLO massless matrix elements. 
The most important outcome of this study is the reduced sensitivity of the spectrum on the choice 
of renormalization and factorization scales. 
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2 Charmonium Production 



Charmonium production in hadronic collisions has been measured both by UA1 |27| and CDF 



0, and has been studied within QCD for quite some time now p8| [2T]. However, all of the 
results available for production at large pt are relative to leading order calculations. Therefore the 
theory of quarkonium production is not as solid as the theory of open heavy quark production, 
and unfortunately a complete NLO calculation presents serious difficulties and lies years ahead of 
us. 

At leading order, large pt charmonium is mostly produced via gluon fusion diagrams such as 
the one shown in fig. |^. Simple power counting in the propagators suggests that the production at 
large p t is strongly suppressed: da /dp 2 ~ l/pt 8 - This suppression can be understood physically 
as a sort of form factor of the charmonium state, form factor which is probed at Q 2 values of the 
order of p t : it is hard to accelerate a J/ip to a large p t without destroying it in the process! 

It was recently realised by Braaten and Yuan [29| that there is a class of higher order processes 
which are not suppressed at large p t as the leading order ones. These are contributions where a 
high p t gluon or charm quark will fragment into charmonium states (Fig. 0). While the probability 
that such a parton fragments into, say, a J/tp is very small, this probability is a number relatively 
constant with p t and the J / ip will carry away a large fraction of the parton momentum. Therefore 
the p t spectrum of charmonium states produced according to this fragmentation mechanism is 
similar to that of the initial hard parton, namely something with an approximate shape da /dp 2 ~ 
l/pt 4 - 

Fragmentation functions for charmonium production from gluon and charm quarks have re- 
cently been calculated and their scaling violation studied |[29|| . They can be convoluted with the 
inclusive p t distributions of gluons and charm to give inclusive p t distributions of charmonium 
states P0f . Some pieces are still missing, such as the fragmentation production of x states. Work 
on these is in progress, and for the time being a final prediction for the spectrum of J/^'s is still 
unavailable. 

A complete fragmentation calculation is available however for the ip(2S), which does not receive 
contributions from x decays. In fig. I show the ip(2S) pt distributions predicted at 1.8 TeV from 
the three processes: direct LO production, gluon fragmentation and charm fragmentation. The 
fragmentation channels become larger than the LO contribution as soon as pt > 7 GeV. In spite 
of this, however, the sum of all these terms is still small compared to the effect of B decays, and 
small compared to the current discrepancy observed between CDF data |7j and the sum of B 
decays and direct production. A similar problem is apparent from the preliminary studies of the 
B fraction in J/ip data [|3l], indicating a large contribution from non-B decays, incompatible with 
current estimates of direct J/ip production. 

It is clear that more has to be learned about charmonium production in hadronic collisions. For 
example, there is a class of processes in between the LO and the fragmentation ones, which behaves 
like l/pt 6 and which probably contributes significatly to the p t region below 10-15 GeV. These 
processes have not been calculated yet. New data and analyses dedicated to separate all possible 
sources of charmonium, for example by separating samples on the basis of the J/ip isolation and 
impact parameter, will be fundamental to help theorists disentangle what today seems a rather 
embarassing issue for QCD. 
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Figure 3: Comparison between different contributions 
to 1p{2S) production. 



Figure 4: Azimuthai correlations between b pairs, as 
a function of different triggers. 
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3 Heavy Quark Jets 



When discussing heavy quark production in hadronic collisions, is has become customary to refer 
to properties of the heavy quark themselves rather than to the properties of the jets in which 
they are embedded. For this reason, available calculations are always formulated in terms of the 
Pt spectrum of bottom mesons and not of their jets. Nevertheless, several measures of current 
interest, say tagging of b-jets produced in association with W bosons, are formulated in terms 
b-jet energies. What makes this distinction interesting is the fact that different mechanisms will 
dominate production of B mesons at a given E t and production of b-flavoured jets of the same E t . 

With a language which is more pictorial than formal, we usually identify three processes for 
heavy quark production: direct production, induced by gg — > QQ scattering, gluon splitting, given 
by a standard QCD gg — > gg collision followed by the branching of a final state gluon into a QQ 
pair, and flavour excitation. The frequently asked question "which is the dominant process for b 
production?" admits several answers, depending on the details of the events we are considering. 

For example, say we are interested in the azimuthal correlations, A(j>. If we consider a sample 
of events triggered by a 20 GeV b, we expect that phase space will suppress in part configurations 
at A(f>=0 relative to A0 = n. In fact configurations with A0=O require the presence of a stiff 
gluon jet recoiling against the bb pair, in order to compensate not only the momentum of the 
triggered b, but also the momentum of the b flying parallel to the b. We expect the overall A<p 
distribution to be peaked at A0 = n. 

If we trigger instead on a 20 GeV b-jet, namely a 20 GeV jet which contains a b quark, 
regardless of its momentum, configurations with A0=O and A<p = n have the same phase space, 
as in both cases it is sufficient the presence of just two 20 GeV jets: in one case the first jet 
is made of the bb pair and the second jet is a gluon; in the other both jets are b jets from the 
direct production channel. The dynamics of the collinear gluon splitting enhances the A0=O 
region, compensating the presence of an additional power of a s , and one expects the resulting 
distribution to have a double peak structure, with peaks at A<j)=0 and tt. 

All of these considerations are qualitative, and will vary with energy and other boundary 
conditions. Nevertheless one can observe their relevance by looking at fig. |], where the two 
different cases are displayed for 20 and 50 GeV jets. The calculation was performed using the 
fully exclusive NLO calculation by MNR which allows to define b-jets using the standard 



cone algorithms applied to inclusive jets in hadronic collisions. Here we define "jet" a set of one 
or two partonsi with AR < 0.7 (AR = V Arj 2 + A0 2 ). More specifically a b-jet is a jet in which 
at least one of the forming partons is a b, while a gluon jet is a jet just made out of one gluon. 
At the NLO, the p t of a gluon jet is by construction equivalent to the p t of the heavy quark pair. 

Using the work of ref . ||13|| , one can calculate the b-jet inclusive E t distribution, and then 
answer at the NLO in QCD another important question, namely "what is the fraction of jets of 
a given E t containing a heavy quark". To do this it is sufficient to compare with the available 
NLO calculations of the jet inclusive E t distributions ||32|| . Fig. |5| plots three NLO quantities: the 



inclusive b (and b) p t distribution (i.e. the NDE spectrum), the gluon jet p t distribution and the 
b-jet inclusive E t distribution. We applied no r\ cuts in these plots. Taking the ratio of the heavy 
quark jet E t distribution with the inclusive jet E t distribution, we obtain the curves shown in 
fig. H for charm and bottom jets. As a denominator, we used a LO calculation for the inclusive 
jet E t spectrum. For AR = 0.7 and fi = p t /2 this was shown to lead to a result numerically equal 



2 At the NLO partonic jets contain at most 2 partons, as there are no more than 3 partons in the final state. 
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Figure 5: Inclusive p t distribution of b quarks, 6-jets 
and gluon jets in b events. 



Figure 6: Fraction of jets containing a b as a function 
of jet E t (solid); ratio between the b Pt and inclusive 
jet Et spectra. 
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to the NLO result. The results are consistent with indipendent estimates of the heavy quark 
content of gluon jets [33|, and are consistent with calculations performed |33| using the shower 



Monte Carlo HERWIG [pq] . The flattening of the fractions at large E t reflects the diminuishing 
percentage of gluonic jets in the inclusive jet sample with larger E t . 

Since, as observed above, measurements based on the b pt or b-jet E t emphasize differently the 
various production mechanisms, a direct measurement of the b-jet E t distribution would provide 
a check of the NLO calculations which is complementary to the one performed using the standard 
NDE spectrum and the B meson p t distribution. 

3.1 Heavy Quark Jets and Photons 

Associated production of photons and charm quarks has been suggested as a probe to study the 
charm density inside the proton [p3fl . The dominant production process is in fact gc — > 7c. The 
similar process gb —>■ ^b would allow the study of the b density. Additional interest in this process 
arises because it is a potential background to the observation of the neutral current decay of a 



4-th generation J 3 = —1/2 quark, b' — > by p6| . Finally, these processes provide a useful control 
sample for the study of heavy quark tagging. 

These processes are theoretically interesting because the density inside the proton of a heavy 
quark is in principle calculable perturbatively ||371| . Neglecting higher order logarithmic corrections, 
which can be resummed using the Altarelli Parisi evolution, the inclusive process pp —>■ 7Q + X 
can be calculated by evaluating the partonic process gg — > QQ-y and integrating over the phase 
space of the Q. This process is dominated by configurations where the quark being integrated over 
is produced at large rapidity and small p t . No divergence will appear, because of the heavy quark 
mass. A consistent definition of the heavy quark density, including thresholds effects, should then 
reproduce the result of this calculation. Comparison against the experimental data is however 
important, to verify that no additional non-perturbative effect is at work. 

I show in fig. |7| the photon p t distribution in events with a central b, calculated using the 
structure function approach (bg — > by, solid line, CTEQ PDF's [|T^]) and the exact matrix elements 
for the gg + qq —> iQQ processes, where higher order corrections to the initial state evolution have 
not been included. The qq channel produces a heavy quark pair via gluon splitting, and cannot be 
accounted for by the structure function calculation. As the plot shows there is perfect agreement 
between the two approaches, at least in the region where gluon splitting is small. This indicates 
that the effects of initial state evolution for the b quark at these values of x and Q 2 are not 
important. Notice that this conclusion is important for the consistency of the NLO evaluation of 
the b cross sections, which is performed without including flavour excitation diagrams [0], already 
accounted for at this order in perturbation theory by the NLO matrix elements. Similar results 
are obtained in the case of charm production, where however there is some larger sensitivity to 
the choice of the charm mass. 

Fig. ^finally contains the distributions of jets of various type produced with photons. Because 
of the difference in charge and partonic density, associated production of b is suppressed by a 
factor equal to 8-10 w.r.t. charm. As a curiosity, notice that because of the suppression of the 
light-quark annihilation channel, it is more likely for a jet produced in association with a photon 
to be a charm jet rather than a gluon jet, at least for trasverse momenta up to 30 GeV. 
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Figure 7: Photon p t distribution in events with a cen- 
tral 6-jet in direct photon events, showing the results 
of the LO structure function approach and of the NLO 
calculation. 



Figure 8: Jet-type composition in 7-jet events. 
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3.2 Associated Production of W and 6-jets 



This is a hot subject in these days, considering its implications on the search for the top quark. 



A detailed study can be found in ref. | 3q| . I will present here some additional remarks, mostly 



matured from the work of and discussions with various CDF collaborators, in particular G. Unal 



and B. Williams 34 



Given the importance of top search, it is important to be able to perform estimates of the 
backgrounds which are as much as possible independent of absolute normalizations provided by a 
theoretical calculation. The most straightforward approach would be therefore to simply assume 
that the fraction of b-jets in the inclusive jet sample is the same as in the Vy+jets sample. The 
first can be measured, and the rate as a function of jet E t and possibly other jet features can be 
then applied to the PF+jet sample (see Contreras presentation at this Workshop). In this fashion 
any theoretical uncertainty on the absolute rate of inclusive VT+multi-jet would drop out, being 
replaced by the observed cross sections. 

A priori there are several reasons to expect this assumption to be too naive. First of all the 
inclusive jet sample has a large contribution from direct production of b jets, contribution which 



is absent in the case of the W sample, where all heavy quarks come from gluon splitting [|3£ 
Secondly, the average fraction of gluon-initiated final state jets is different in the two processes. 
This difference will affect the rate of heavy quarks produced via gluon splitting. Both these reasons 
point to a smaller fraction of heavy quark jets in the W sample than in the inclusive jet sample, 
indicating that at worst the above background assumption was conservative. Finally, even gluon 
jets of the same E t in the two samples will have a different probability to produce QQ pairs, 
because of the intrinsically larger Q 2 scale of the collision producing a massive W . 

In spite of these differences, we expect the calculation of the relative fractions of heavy quark 
jets to be less affected by theoretical uncertainties than the absolute rates. The predictions can 
be tested on the inclusive sample of W+ one 6-jet + X, which is not significantly contaminated 
by top. Notice that for this calculation it is fundamental to keep into account the effects of the b 
mass, in order to avoid divergencies associated with the integration over the second b phase space 
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Without getting into the details of the calculation, performed using a mixture of the matrix 
element results of ref. |38| and HERWIG | 35fl , I show an indicative result in fig. |9|, which gives the 



fraction of iy+multijet events containing a 6-jet above 20 GeV, as a function of the multiplicity 
of jets above a given E t threshold. It appears that the probability per jet of being a 6-jet is of 
the order of 1%, almost independently of the jet multiplicity. For comparison, I show in fig. ITD 



the same distribution obtained for generic multi-jet events. As expected, the fractions here are 
slightly larger, indicating that indeed the naive assumption that jets in the two samples have the 
same 6-quark content is conservative. 

Acknowledgements: I wish to thank my collaborators on heavy quark physics, namely U. Baur, 
E. Braaten, M. Doncheski, S. Fleming, S. Frixione, S. Keller, P. Nason and G. Ridolfi. I also wish 
to thank my collaborators at CDF for the constant stimulus provided by their terrific work. 
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Figure 9: Fraction of VF+multijet events containing a 
fe-jet with E t >20 GeV, as a function of the number of 
jets above a given threshold. 



Figure 10: Same for inclusive multi-jet events. All jets 
have \r)\ < 2. 



12 



[1] P. Nason, 'Heavy Quark Production', in "Heavy Flavours", Advanced Series on Directions in High 

Energy Physics, ed. M Lindner, AJ Buras. Singapore: World Scientific. 
[2] C. Albajar et al., UAl Coll., Phys. Lett. 213B (1988), 405. 

[3] P. Nason, S. Dawson and R. K. Ellis, Nucl. Phys. B303 (1988), 607; Nucl. Phys. B327 (1988), 49 . 
[4] C. Albajar et al., UAl Coll., Phys. Lett. 256B (1991), 121. 

[5] W. Beenakker, H. Kuijf, W.L. van Neerven and J. Smith, Phys. Rev. D40 (1989), 54; W. Beenakker, 

W.L. van Neerven, R. Meng, G.A. Schuler and J. Smith, Nucl. Phys. B351 (1991), 507. 
[6] F. Abe et al., CDF Coll., Phys. Rev. Lett. 68 (1992), 3403; Phys. Rev. Lett. 71 (1993), 500; 

FERMILAB-PUB-93-106-E; FERMILAB-PUB-93-145-E. 
[7] F. Abe et al., CDF Coll., Phys. Rev. Lett. 69 (1992), 3704. 
[8] J.C. Collins and R.K. Ellis, Nucl. Phys. B360 (1991), 3; 

S. Catani, M. Ciafaloni and F. Hautmann, Nucl. Phys. B366 (1991), 135; 

E.M. Levin, M.G. Ryskin and Yu.M. Shabelsky, Phys. Lett. 260B (1991), 429. 
[9] E.L. Berger, R. Meng and W.K. Tung, Phys. Rev. D46 (1992), 1895; 

E.L. Berger, R. Meng and J. Qiu, ANL-HEP-CP-92-79. 
[10] M.L. Mangano, Z. Phys. C58 (1993), 651. 

[11] A. Martin, R. Roberts and J. Stirling, Phys. Lett. 306B (1993), 145. 
[12] J. Botts et al, Phys. Lett. 304B (1993), 159. 

[13] M. Mangano, P. Nason and G. Ridolfi, JVucI. Phys. B373 (1992), 295. 

[14] A. Geiser, UAl, this Workshop and Proceedings of the XXVII Rencontres de Moriond, Les Arcs, 
Perturbative QCD and Hadronic Interactions, ed. J. Tran Thanh Van, (1992), p. 159; PhD thesis. 
RWTH, Aachen (1992); 

C. Albajar et al., UAl, CERN-PPE/93-153, to appear on Z.Phys. C. 
[15] F. Abe et al., CDF Coll., FERMILAB-CONF-93-199-E, submitted to the 16th Intern. Symp. on 

Lept. Phot. Inter., Ithaca NY, Aug. 1993; T. Huffman, CDF Coll., FERMILAB-CONF-93-274-E. 
[16] E. Berger and R. Meng, CERN-TH.6950/93 

[17] E. Laenen, J. Smith and W.L. van Neerven, Nucl. Phys. B369 (1992), 543. 
[18] F.A. Berends, J.B. Tausk and W.T. Giele, Phys. Rev. D47 (1993), 2746; 

E. Laenen, J. Smith and W.L. van Neerven, FERMILAB-Pub-93/270-T (1993). 
[19] J. Skarha, A. Holscher and D. Crane, CDF, presentations at this Workshop. 
[20] A. Maciel and S. Igarashi, DO Coll., presented at this Workshop. 
[21] E.W.N. Glover, A.D. Martin and W.J. Stirling, Z. Phys. C38 (1988), 473; 

E.W.N. Glover, F. Halzen and A.D. Martin, Phys. Lett. 185B (1987), 441. 
[22] See M. Contreras (CDF), B. Denby (CDF) and H. Greenlee (DO) presentations at this Workshop. 
[23] R.S. Fletcher, F. Halzen and E. Zas, Phys. Lett. 221B (1989), 403; 

U. Baur et al, Florida State FSU-HEP-930816, to appear in Phys. Lett. B. 
[24] M. Mangano, in preparation. 

[25] V.A. Khoze, J. Ohnemus and W.J. Stirling, Durham preprint, DTP/93/52; 
V.A. Khoze, L.H. Orr and W.J. Stirling, Nucl. Phys. B378 (1992), 413; 

for a review of this and related topics in heavy quark QCD, see: V.A. Khoze, Durham preprint 
DTP/93/78. 

[26] M. Cacciari and M. Greco, Univ. of Pavia FNT/T-93/43, [hep-ph/93 11260 . 
[27] C. Albajar et al., UAl Coll., Phys. Lett. 256B (1991), 112. 
[28] E.L. Berger and D. Jones, Phys. Rev. D23 (1981), 1521; 

R. Baier and R. Riickl, Z. Phys. C19 (1983), 251; 

B. Humpert, Phys. Lett. 184B (1987), 105; 

R. Gastmans, W. Troost and T.T. Wu, Nucl. Phys. B291 (1987), 731. 
[29] E. Braaten and T.C. Yuan, Phys. Rev. Lett. 71 (1993), 1673. 



13 



[30] M. Doncheski, S. Fleming and M. Mangano, FNAL-CONF-93/348-T. 
[31] F. Abe et al., CDF, FERMIL AB-PUB-93- 1 58-E. 

[32] F. Aversa et al., Phys. Lett. 210B (1988), 225; Z. Phys. C49 (459), 1991; 

S. Ellis, Z. Kunszt, D. Soper, Phys. Rev. Lett. 62 (1989), 2188; Phys. Rev. Lett. 64 (1990), 2121. 
[33] M.L. Mangano and P. Nason, Phys. Lett. 285B (1992), 160. 
[34] G. Unal and H. Williams, private communication. 
[35] G. Marchesini and B. Webber, Nucl. Phys. B310 (1988), 461. 
[36] W.S. Hou and R.G. Stuart, Phys. Rev. Lett. 62 (1989), 617. 
[37] J.C. Collins and W.K. Tung, Nucl. Phys. B278 (1986), 934. 
[38] M.L. Mangano, Nucl. Phys. B405 (1993), 536; 

Z. Kunszt, Nucl. Phys. B247 (1984), 339. 



14 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/93 1 1378vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/93 1 1378vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/93 1 1378vl 



